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Abstract 

We present an intuitive physical picture for fractures in nacre-type stratified materials via scaling 
arguments: strain distributions around a fracture are rather different depending on directions and 
size of fractures. We thus observe that viscoelastic effects are important for parallel fractures. This 
effect are taken into account via the simplest viscoelastic model for weakly cross-linked polymer. 
Within a certain limit, we find a trumpet crack shape similar to that in certain polymers and make 
predictions for the fracture energy of certain stratified materials in this mode. 
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I. INTRODUCTION 



Strong structures created in nature through a long history of natural selection sometimes 
take advantage of composite structures such as laminar structures (tooth, timber etc.) where 
soft and hard layers are intertwined. Biomimetics aims at exploiting such ideas and some 
industrial materials (from blockboard to high-tech tires) have utilized the concept of stacking 
lamellae to enhance their toughness [I], 0, Q |], In a series of papers || [7|, |8|, we have 
studied another example of such soft-hard laminar structure on nanoscales, nacre ||, |IU|. [IT], 
to present physical pictures for this substance. 



perpendicular fracture 




parallel fracture 



FIG. 1: Narce-type structure of materials: hard layers (the elastic modulus Eh) has a typical 
thickness dh, which are glued together by soft layers (the elastic modulus E s ) of a typical thickness 
d s . An important conditions for these parameters are Eq = E s /Eh,£d = d s /dh -C 1 with e = 
£ o/ e d *C 1. The cracks in the y — z plane and the x — z plane are called a perpendicular and 
parallel fractures, respectively (the z-axis is perpendicular to this page). 

The structure and notations are illustrated in Fig. [1]. A typical thickness d s of the soft 
layer made from a protein, conchiolin, is dozens of nanometer, which is much smaller than 
a micrometric thickness d s of the hard layer composed of aragonite. An elastic modulus E s 
of the soft layer is also much smaller than Eh of the hard layer where the latter is roughly 
given by 50 GPa. We define two small quantities: 

e d = d s /d h (1) 
e = E s /E h (2) 

An important property of nacre is that e defined by 

£ = e /e d (3) 
2 



is very small. We will work in the limit of Eo,Ed,s <C 1, and also in the continuum limit 
where relevant length scales are much larger than the layer spacing d: 



d = d s + d h . (4) 

We emphasize here that under these conditions our theory is applicable to a general layered 
structure other than nacre. Physically, Ed has to be small to make a relatively ridged material 
(as a whole), while e (and thus Eq) has to be small to make the stress concentration around 
tip small to enhance the toughness] these conditions seem to have been indeed selected in 
nacre by nature. 

We have been studied both perpendicular and parallel fractures (see Fig. P within the 
linear elastic fracture mechanics (LEFM |13| ) under the plain strain condition (the sample 



is thick in the ^-direction) |6], [7], [8] . One of our aims here is to include a linear viscoelastic 
effect into this system. 



II. NOVEL VIEWS ON PREVIOUS RESULTS 

The proposed elastic energy per unit volume of the stratified system is given by @ 

/ = E h(? xx + E e 2 yy + E e 2 xy + E e xx e yy (5) 

where 

E = eE h , (6) 

We have introduced the strain field 

for the deformation field Ui where (xi,x 2 ,x 3 ) = (x,y,z). This energy reflects, via = 
df/deij, intuitively self-evident properties that (1) this system is hard for a tensile stress 
in the x-direction (a xx ~ Ehe xx ) due to hard layers, (2) it is soft for a tensile stress in the 
y-direction (a yy ~ E e yy ) due to soft layers, and (3) it is soft for a shear stress [a xy ~ E Q e xy ) 
again due to soft layers. Here, we notice that the effectiveness of soft layers is diminished 
(E a — > Eq) by the thin thickness of the soft layers. 

For parallel fractures, we have shown that a usually negligible bending term becomes 
important ||, which shall be confirmed below from a different viewpoint. The bending term 
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is given by 



f B 



K B 




(8) 



where 



K B 



(9) 



Here we have introduced a length scale, 



I = dj y/e, 



(10) 



which is much longer than the layer period d. This bending term becomes important also 
in certain liquid crystals |14], However, there is a significant difference from the present 
situation; in liquid crystals I is replaced by a much smaller length of the order of atomic 
scales, which results in a strain distribution quite different from nacre (see below). 

A. Scaling structure of elastic energy 

An important difference from the conventional isotropic elastic theory is that, in an 
anisotropic system such as nacre, there exist two distinct length scales X and Y . On the 
contrary, in isotropic systems at equilibrium when the elastic body is deformed over a range 
of R in the x-direction, the deformation relaxes out at the same distance R also in the 
y-direction; this is because the deformation field of an isotropic system satisfies an elliptic 
differential equation with a unique length scale (under the plain strain condition). In other 
words, the deformation in nacre does not satisfy the Laplace equation. Note here that, in 
the context of a fracture problem, X corresponds to a length of parallel fracture while Y to 
that of perpendicular fracture (see Fig. || below). 

Thus, a local energy (per unit volume) of nacre is dimensionally expressed as 



At equilibrium (statics), / is locally minimized for both variables u x and u y ; we have two 
equations for u x and u y , i.e. df /du x = and Of /du y = 0. Seeking the relevant solution by 
requiring (u x , u y ) ^ (0, 0), we have an equation for X and Y . This equation has two solution 
specifying a relation between X and Y; one corresponds to perpendicular fractures while 




(11) 
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the other to parallel fractures. We should note that here and hereafter a limitation of this 
type of dimensional expressions: relative signs (or phases of complex numbers) sometimes 
can not be specified. 



1. Parallel fractures 

The solution for parallel fractures announces a relation, X 2 ~ Y 2 (l + l 2 /X 2 ), which 
implies two regimes for this problem: X ^> I and X -C Z. 
For large fractures (X 3> I), we have an isotropic relation, 

X ~ Y. (12) 

On the other hand, with an aide of a relation df/du x = together with Eq. (|T2"D, we obtain 
an anisotropic deformation field, 

u x ~ eu y ] (13) 



we find u x <C u y , which is appropriate for parallel fractures. Using Eqs. (|i2[ ) and fll3D , we 
can estimate magnitudes of each term in Eq. ([0]). Keeping terms only at the leading order 
in e, we arrive at 

/-MI) 2 < 14 > 

where we denote scales of X (~ Y") and w y as i? and w, respectively. This corresponds to a 
more precise form / ~ E (^jffj + E (^jr^J shown in ]7[. 
For small fractures (X <C i) , we have instead 

x 2 ~ yz, (15) 

which implies X ^> F. The deformation field satisfies 

Ux ^ (16) 

Here, we note an important condition for a continuum theory: X,Y > d. This condition, 
together with Eq. (|10D, results in 

sl 2 /X 2 ,el 2 /Y 2 <1 (17) 

which implies el/X,el/Y < ^fe <ti 1 (note also that this does not necessarily imply I < X,Y; 
e.g. X > y/el satisfies both I > X and el/X < 1). Thus, we again find u x u y . In this 



way, we arrive at 

2 



f*Eo(£) *K B (£) (18) 
which corresponds to / ~ E Q (^fj + Kb (tjJt) shown in ||. 



2. Perpendicular fractures 

The solution for perpendicular fractures announces an anisotropic relations, 

Y ~ ^feX (19) 
u y ~ y/eu x (20) 



for which u x ^> u y . These relations lead to an energy, 



Here, u denotes not u y but u x . This energy corresponds to / ~ E h [^-) 2 + E (^^j shown 



in 



B. Scaling views on fracture mechanics 

We consider a fracture in nacre-type materials. Strain distribution for three types of 
fractures are schematically given in Fig. |2|. We should note that these figures are at the 
scaling level and thus they are rather rough images; some examples of comparison with 
exact results are given in Appendix. We will develop below scaling arguments in order 
to reproduce our previous results but only dimensionally. Similar arguments in a different 



context was first presented in [14 



1. Parallel fractures 

For large cracks (X ^> /), the strain distribution is isotropic (Fig. 0a) and, thus, the 
potential energy per unit crack-front length (in z-direction) is dimensionally given by 

F ~ ^o(^) 2 R 2 - cruR + 2G\\R (22) 
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(a) Large parallel fracture 




(b) Small parallel fracture 




(c) Perpendicular fracture 




FIG. 2: Physical images of strain distribution accompanying a fracture. Dark areas indicate the 
region where the strain field is significant, (a) Large parallel fracture. There is only one length 
scale R (S> I) for the distribution, where the crack length is also of the order of R. (b) Small 
parallel fracture. The strain distribution becomes rather anisotropic: the length scale X (<C I) for 
rc-direction is larger than Y, where the crack size is of the order of X. (c) Perpendicular fracture. 
An anisotropic strain distribution with I>7, where the crack size is of the order of Y. 

where the first two terms describes the elastic potential energy and the last term an energy 
loss due to a creation of new surfaces in this system: G» is the fracture energy (energy 
required to create a unit area). We minimize this energy with respect to u to find a Hooke's 
law 

a ~ Equ/R (23) 
which leads to an optimized energy value, 

F m ~ -cr 2 R 2 /E + 2G\\R, (24) 

which is quadratic in R; the maximum is given at R = R*: 

a ~ y/G\\E /R*. (25) 

When R < R*, F m decreases with decrease in R; a fracture with the size smaller than R* 
tends to close to lower the energy. On the contrary, when R > R*, F m decreases with 
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increase in R; a fracture with the size larger than R* tends to expand. Thus, Eq. fl25|) 
corresponds to a critical failure stress. From Eqs. (|23| ) and (f25|), we have 




Eqs. fl2"oD and ( |2"6| ) give scaling structures for the stress and deformation fields obtained 
more precisely in J7|: the stress scales as R~ x l 2 and the deformation is parabolic ^~ VR^j 
as in the conventional LEFM. They also give a scaling relation 



au ~ G 



(27) 



which announces that the product au (at critical state, R ~ R*) gives the fracture energy 
(see Discussion). The fracture energy G\\ might be associated with the separation energy 
of the soft layer and we previously concluded that there was no significant enhancement in 
this mode. Due to the ensuing analysis we can afford a more precise understanding on this. 

In conventional LEFM, the stress intensity factor K scales as \JQEq in the present context 
where Q is the energy release rate and its critical value is the fracture energy (G\\). Eq. fl25|) 
suggests that K scales not as \/QE but as its critical value ^JG\\Eq] this is logical because 
in Eqs. (|2^|)-(f27|) we are always at a critical of failure (R ~ R*). 

For small cracks (X <C I), the strain distribution is anisotropic (Fig. ^|b). We note here 
that this distribution is quite different from certain cases of liquid crystals |TJ|]; there strain 
is distributed more widely in y-direction than in x-direction. The potential of our system is 
given by 



which corresponds to a yy ~ E du y /dy. Note here that the dominant component of stress 
field tensor is a yy (e.g. o xy ~ E du y /dx ~ E (u/X) <^ o yy due to Eq. (p~5|) ) . The energy 
optimized for u is here given by 




(28) 



We minimize this energy with respect to u to find a Hooke's law 



E u/Y 



(29) 



F, 



m — 



1 1 K B + 2G\\X. 



(30) 
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F m takes its maximum value at X ~ X* ~ \/Y*l, where 



//en , /Guy* , N 

Eqs. ( |3~1~D and give scaling structures for the stress and deformation fields as well 
as a scaling relation au ~ GV We emphasize here that the stress thus obtained (which is 
proportional to strain) is indeed anisotropic; <r~l/X~l/ \fY~l. The stress and deformation 
fields obtained in || are also consistent with these scaling structures; we can check that, 
dimensionally, a(x,y) and u(x,y) derived in || indeed reduce to Eqs. (|31|) and (|32| ) by 
setting x 2 ~ (because we are always in the regime specified by Eq. fll5|) at the scaling 
level) . 

2. Perpendicular fractures 

The strain distribution is anisotropic (Fig. 0c) and, thus, the potential is given by 

F~E h [^j 2 XY - auY + 2G ± Y } (33) 

where G± is the fracture energy. We minimize this energy with respect to u to find a Hooke's 
law 

a ~ E h u/X (34) 

which corresponds to o xx ~ Ehdu x /dx (the dominant stress tensor component in this case 
is u xx ). The energy optimized for u is here given by 

a 2 Y 2 

F m ~-—— + 2G ± Y. (35) 
F m takes its maximum value at Y ~ Y* ~ y^X*, where 
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(36) 



n^J^(r*/Vi), (37) 

Scaling structures for the stress and deformation fields in Eqs. (|36D and fl37|) are in ac- 
cord with results in (the stress or strain thus obtained is indeed anisotropic). In this 



case, we can estimate the fracture energy au ~ G± via the maximum stress for a con- 
tinuum stress (ctf — y/ yfeG±Eh/dJ balancing with the yield stress of the pure aragonite 
(ays — \/ Eh^h/ a-hj '■ G± ~ X±jh where Aj_ = dj (yfeah) 1; the fracture energy is en- 
hanced from that of the pure aragonite 7^ because ah (-C d) is a microscopic size of defects 
in a sample. 

III. INCLUSION OF VISCOELASTIC EFFECT 

From the above scaling arguments, we have seen that Vby IS Ct dominant component of 
the deformation field for parallel fractures while u x is important for perpendicular fractures. 
Since u y essentially corresponds to the deformation in soft layers while u x to that in hard 
layers, viscoelastic effects associated with extremely soft layers becomes important only for 
parallel fractures; we shall include these effects into our theory. 

As a model for viscoelasticity, we employ the simplest model (also utilized in [[l6|j ): 

fi(u) = /2 + (/Ax. - Ho) ; l —r (38) 

1 + IUOT 

This is the minimal model to include the following essential properties of weakly cross- 
linked network: (1) its slow motions are governed by a weak-modulus fiQ associated with 
weak cross-links, (2) its fast motions involve a strong-modulus yUoo fio) originating from 
entanglements. 

a=ai+a 3 



a 3 =E 3 e 3 



e=ei+e 2 =e 3 
CTi=a 2 




FIG. 3: Zener model with two springs and a dashpot. 



This model is a special case of the Zener model (see Fig. ^|), |T7| , |18|1 for which stress a 
and strain e satisfy the relation, 

^+> = £ 3 e + ?A^e (39) 
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with identification, E 3 = /i , E 2 = fioo — Ho, and r = r]/E 2 . 

The previous purely elastic treatment should correspond to the zero frequency limit or 
the static limit of a viscoelastic model. For parallel fractures, (y, y) component of stress 
(a) is important, for which we write a = a s ~ E s e s where a s and e s denote stress and 
strain of soft layers. Thus, it is natural to require that soft layers obey the trumpet relation 
a s + T& s = fi e s + /iooTe s with /i = E s , because then the static limit reduces to the desired 
form. Noting the relations between coarse-grained stress and strain of nacre (a and e) and 
more local stress and strain of soft layers (cr, and e s ), i.e. o = o s and e s d s = ed, we deduce 
a viscoelastic model for nacre: 

a + rcr = E e + E^re (40) 

where 

= XE = XE s /e d (41) 
r = (42) 

A = /Uoo//io (43) 

We note again that fio = E s is associated with weak cross-links in soft layers, while fioc, with 
entanglement in soft layers. 

Putting a = <j m e iujt and e = e m e tujt , with the definition of the viscoelastic modulus, 
c"m = Mci;)e m , we obtain Eq. ( |38| ) with and /i replaced by E^ and E , respectively. 
This model can be divided into three regimes: 

I. ujt 1/A: /i(co>) ~ E (soft solid) 

II. 1/A < cur < 1: ix{uj) ~ iuJT) (liquid) (44) 

III. 1 ujt: fi(uo) ~ Eoo (hard solid) 

IV. CRACK SHAPE AND FRACTURE ENERGY 

In the following we consider only large cracks: / ^> X. The opposite limit is examined in 
a separate article. [19|] In the case of such large cracks, we are always in the isotropic regime 
where X ~ Y{= R) (see also Fig. 0a). We imagine a crack propagating at a constant speed 

V. For simplicity, we treat only certain large velocities: / Vt; although it is very difficult 
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to estimate r for the protein, this corresponds to r > 10~ 2 s for / ~ 100 /im and V ~ 1 
cm/s. 

An important observation here is the scaling identification of a distance r from the tip 
and the frequency uj via a stationary crack-propagation speed V: 

r ~ V/w; (45) 

small distances correspond to high frequencies while long distances to low frequencies — the 
farther away from the tip, the more time for relaxation. Corresponding to Eq. (fH|), the 
fracture is spatially divided into three regions: 

I. AVr <C r: soft solid 

II. yr«r« AVr: liquid (46) 

III. d < r < Vr: hard solid 

Another important observation is that, in certain model of viscoelasticity, the scaling 
relation a ~ r^ 1 / 2 (where r is a distance from the fracture tip) still holds as in a purely 
elastic model. This point shall be explained in next section. 

First we consider a large crack size (L ^> XVt); then, all the three regions I-III are 
present. The soft-solid region I corresponds to low frequencies, and thus, to the static limit; 
in this region (r ^> XVt ^> I), Eqs. (|25|)-(|2"T|) holds. In the liquid zone II, the stress field 
scales as o ~ ur/e ~ r/Vu/r 2 and at the same time it should scale as a ~ r -1 / 2 from the 
above second observation. Thus, the strain should scale as u ~ r 3 / 2 and the product as 
ou ~ r. The coefficients can be determined by matching at r ~ XVt, i.e. the boundary 
between I and II: for example, we have 

T 

ou — G\\ (liquid zone). (47) 

XV T 

When we reach the hard solid region (r ~ Vr), we find 

o"M~G||/A (hard-solid zone) (48) 

When we are in the solid region, as in the same manner (due to the simplifying condition 
/ Vr) with the derivation of Eqs. (p5|)-(|27|), we can show 



a ~ ^GoE^/R* (49) 



u ~ y/GoR'/Eoo (50) 
ou ~ G (51) 
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Here, we emphasize that Go is associated with the hard solid appearing near the tip and is 
different from Gy (At the place very close to the tip (r <C I), Eqs. (f4l^-(|5T|) are no longer 
valid and they merge to the lenticular expressions ||). From Eqs. (|i8D and (|51"|), we find 



Gii ~ AG . 



(52) 



The overall separation energy G\\ is enhanced from Go associated with the hard solid. The 
crack shape resulting from this analysis is just like a trumpet (see Fig. |j), as has been 
suggested by the name of the model. The similar form has been predicted |16|] and observed 
20fl also in systems (weakly or non cross-linked polymers) rather different from the present 



anisotropic system. 




FIG. 4: Viscoelastic trumpet. The dependence of deformation field u on a distance from the tip R 
is indicated in the figure. Very close to the tip (I <C i?), the shape merges to a lenticular form M. 



We complete our arguments by considering smaller fractures. When (/ <)Vt < L < XVt, 
only the hard-solid and liquid regions are present; the soft solid has yet to develop. In 
this situation, the fracture energy is given by Eq. (|47|) at r = L: G\\(V) = GqL/(Vt); 



the toughness decreases with velocity. When (I <)L < Vr, only the hard-solid region is 
developed and the fracture energy is given by Go- Thus, as a function of V, the fracture 
energy starts from a larger plateau value AGo, and then decreases to reach a smaller plateau 
value Go: 

' AG for d/r < V < L/ (Ar) 

G L/(Vt) for Lj (Ar) < V < L/r (53) 
G for L/t < V 

Note here that in a continuum theory we only consider the region Vt > d. This behavior 



Gii GO 
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can be confirmed more precisely from a general formula [16| : 



G\\(V) 



~ E. 



1 



duo 

— 1m 

UJ 



1 



(54) 



Go 



DO 



which can be analytically calculated for the trumpet model (just as in the previously known 
isotropic case [PJ). 

V. DISCUSSION 

In |Tj| the following observations are used: (Al) for a certain range of V, a (dimensional) 
dynamical equation pDv/Dt = — Vct for viscoelastic materials can be replaced by Va = 0, 
(A2) accordingly, in certain model of viscoelasticity, the scaling relation a ~ r -1 / 2 still holds 
as in a purely elastic model, (A3) similar to the usual LEFM, we can expect a relation, 
Q ~ K 2 j 'Eoo for certain viscoelastic models (where K and Q the stress intensity factor and 
the energy release rate associated with the tip), and (A4) we can find scaling structure of 
fracture energy from the product o~u, calculated at a large r. (Our arguments are, however, 
slightly different from the original one in that we did not explicitly used (A3).) We present 
possible interpretations for these observations to complement the original arguments. 

Before showing our interpretations, we assume that the rheological relation is linear in a 
sense that the stress and strain are governed by 



In terms of Fourier transform (in time), this reduces to a simple proportional relation: 



where fJ>(u) = iuE^; (J,(u) has the dimension of a elastic modulus E(t). Here, we have used 
the notation O u = J dtO(t)e- iuJt . We shall use the scaling identification Eq. (|45| ) in the 
following arguments. 

A. Interpretation for V<7 = 

pDv/Dt = — V<7 can be dimensionally reexpressed in terms of Fourier components (in 
space and time) of stress and deformation as pui 2 Uk,u — kok^- This equation reduces to 




(55) 




(56) 
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V 2 Uk,uj — c 2 Uk,uj for a linear rheology (<Tk,w — Ekuk tLtJ ) and for a crack propagation speed V 
(u ~ Vk, or Eq. PS])). Here, the sound velocity c is defined as c = y/Wfp. Thus, if V <C c, 
we can neglect pDv/Dt against V<7. For ~ 10 5 N/m 2 and p ~ 10 3 kg/m 3 , c ~ 10 m/s, 
which is sometimes high compared with typical crack propagation speeds, say, of cm/s. 



B. Interpretation for a ~ r x l 2 

In the plane strain condition, there are only three independent stress components: a xx , 
Uyy and a xy . For V <C c, we can start from Vcr = 0, which gives only two independent 
equations: + = and + = 0. The third equation is a compatibility 
equation, which directly follows from the definition (Eq. (|7])): + = 2-^4^. Again, 
for a linear rheology, we have a simple proportional relation of the type = E^mCm or 
eij = Dijkidki in Fourier space (in time only). Then, the three differential equations for 
Fourier component (in time) of <Jij with respect to variables x and y have the same spacial 
scaling structure as in a purely elastic model where a(t) = E(t)e(t). Thus, the scaling 
structure, a ~ r -1 / 2 should be unaltered even in a linear viscoelastic model. 



C. Interpretation for Q ~ K 2 /E^ and scaling formula for fracture toughness G ~ au 

To establish Q ~ K 2 ^^ in the viscoelastic model and to understand why it can not 
be replaced by Q ~ K 2 /E , we start from the text book derivation [13|]; we imagine that a 



crack is closed to a smaller size from the tip by applying a necessary work, say from a initial 
size a + Aa to a reduced size a (see Fig. |5|). 

The energy release rate Q (per unit length of crack front) is defined by 

G= Urn (57) 

where the work of crack closure (per unit crack-front length) can be expressed as 

px=Aa p rO 

AU ~ / 2dx dy (-a) de (58) 



J x=0 J Je 

Here, the factor 2 counts both the upper and the lower sides of the crack. 

This quantity depends on a history of closure process, for example, via Eq. (|55|). However, 
we will finally consider the limit Aa — > as in Eq. (|57D, while the relevant frequency is 
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y 



(a) initial crack 




(b) reduced crack 

Closure stresses 
at the initial time 




FIG. 5: Closure of a crack. 

given by V/Aa (—>■ oo): very high frequency is important. In such a case, the trumpet model 
behaves as a solid with a modulus E^ (a ~ E^e). Namely, 



AU 



J dxj dy^/E^. (59) 



As shown before, even in the trumpet model, we expect a ~ K/\Jr, if we denote the 
coefficient (the stress intensity factor) by K. Thus, Eq. (|59D leads to AU ~ AaK 2 / E^; we 
indeed recover Q ~ K 2 /E oa . At the critical of failure, this corresponds to 

G ~ i^V^oo (60) 

T/ie critical value of the energy release rate, which is the conventional fracture energy, is 
associated here with the hard solid region near the tip. 

We emphasize here that the overall fracture energy G\\ should scale as liniAa^oo which 
is different from the critical value of the energy release rate. Then, in contrast with the limit 
Aa — ► 0, very low frequency becomes important; we have to use a ~ E e in this context to 
find 

G|| ~ K 2 /E . (61) 



The stress intensity factor K in Eqs. (|60|) and (^T|) are actually different; these factors 
are an unknown coefficient for a known singularity (~ y/r) in the above arguments. In fact, 
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from Eqs. (|25|) and (|49|) , they are respectively given by ^/GqE^ and ^G\\E . In other 
words, Eqs. (^) and (|6lD can be directly reproduced from Eqs. (|25| ) and (f49[). 

Noting AU ~ J (ira^i^a;) for linear models (cr ~ e), we have an expression of the 
overall fracture energy Gn for a macroscopic Aa: 



1 f Aa 

Gll ~Aa~J a ( x ) u ( x ) dx ( 62 ) 

If the product a(x)u(x) reaches a plateau value at a macroscopic distance xo, we can expect 
a relation: 

G\\ ~ (x(x )w(:ro) (63) 

where the left-hand side is evaluated at a macroscopic distance Xq from the tip. This explains 
why we can use this product to estimate the fracture toughness, as announced in Sec. II. B 
(see, for example, Eq. (^7|)). 

VI. CONCLUSION 

In this paper, we present physical pictures for fractures in nacre-type materials via scaling 
arguments, which is complementary to our previous more detailed analysis: strain distri- 
butions around a fracture are significantly different among (small and large) parallel and 
perpendicular fractures. From pictures thus obtained, we see that viscoelastic effects may 
play a role in parallel fractures and the effects are taken into account via the simplest model. 
We limit ourselves to the case of large cracks (X ^> I) where crack-propagation speeds are 
not too slow (I <C Vr); beyond this limit, crack shape and behavior of G(V) seems to be 
quite different, as discussed in a separate paper f]l|. Within this limit, we found that the 



crack shape takes the trumpet form as in weakly (or non-) cross-linked polymers |20] al- 
though they are different in the region very close to the tip (this difference might be difficult 
to detect). The overall fracture energy G\\ is found to be enhanced from separation energy 
Go associated with the hard solid region developed near the tip. It is emphasized that, here, 
the critical value of the energy release rate corresponds not to the overall fracture energy 
but a separation energy associated with the hard solid. 
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APPENDIX: ROUGH IMAGES AT THE SCALING LEVEL AND EXACT DIS- 
TRIBUTIONS 



In order to clarify the meaning of physical pictures illustrated in Fig. |2|, we give examples 
on comparison of such pictures with an exact solution available for perpendicular cracks 
(Figs. |6| and [?p. We employ the following analytical expression J7J: 



a 



Re 



«/£ - 1 



d/ 2 



(A.l) 



with z = x + iy/\/e where e = 1 and e — 0.1 for Figs. || and ^ respectively (L = 1). This 
expression have separate length scales for x- and y-direction: L and Lj 'y/e, respectively. Fig. 
^presents an example of the conventional isotropic sample (e = 1) while Fig. [F|an example 
of a perpendicular crack in a nacre- type material (e = 0.1). 




4 -0.5 0.5 1 

FIG. 6: Left: contour plot of strain distribution for a crack (size L) in an isotropic elastic sample 
(the dimension 2L x 2L) based on an analytical solution (e = 1). Here, L is the unit length. Right: 
corresponding physical image. The distribution is isotropic in a sense that there is only one length 
scale L. 
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-1 -0.5 0.5 1 



FIG. 7: Left: contour plot of strain distribution for a perpendicular crack (size L) in an nacre-type 
elastic sample (the dimension 2L x 2L) based on an analytical solution obtained in |7j (e = 0.1). 
Here again, L is the unit length. Right: corresponding physical image. The distribution is aisotropic 
in a sense that there is two separate length scales where the length scale for x-direction is L and 
that for y-direction is L/y/e. 

In most experimental cases the crack size can be much smaller than the sample dimension. 
In such cases, detail analytical behaviors will change but physical pictures at the scaling level 
stays intact. 
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